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Status of heavy quark physics from the lattice 

Nicolas Garron a 

a DESY, Platanenallee 6, 15738 Zeuthen, Germany 

In this short review, I present a summary of various methods used to simulate heavy quarks on the lattice. I 
mainly focus on effectives theories, and give some physical results. 



1. Introduction 

Lattice simulations provide a powerful tool to 
study QCD beyond perturbation theory. In par- 
ticular, for quark masses around or below that of 
the strange quark, today's lattices offer the pos- 
sibility of very accurate and direct calculations 
of hadronic observables. The recent feasibility of 
simulating dynamical fermions (i.e. without ne- 
glecting the effects of internal quark loops) pro- 
vide a significant reduction (and a better con- 
trol) of the systematic errors. However, simu- 
lating heavy quarks on a lattice is much more 
challenging than simulating the strange quark. 
This is rather unfortunate, because precise calcu- 
lations in the heavy-quarks sector are needed to 
constrain the standard model (see, e.g. [l|2|3|4j ). 
The problem of heavy quarks on the lattice is a 
long-standing one, and various approaches have 
been proposed already more than 20 years ago. 
Of course these ideas have grown up, and together 
with the improvement of the lattice techniques, 
they have become more and more predictive. In 
addition to that, other very interesting strategies 
have been proposed more recently. Some of them 
have already given physical results, and others 
are not in that stage yet. The scope of this talk 
is to present the various methods (and their lim- 
itations) that exist to deal with heavy quarks on 
the lattice. I present also a (limited) selection 
of results which have been obtained by different 
lattice groups. I would like to apologize to all 
colleagues that I cannot cite, due to space limita- 
tions. For an overview of the most recent lattice 
results see [5]. 

This paper is organized as follow: In Sect. 
I review the problem associated with simulating 
heavy quarks on a lattice. In Sect. 02 I present an 



overview of the different approaches. In the last 
section, I report some recent physical results. 

2. Heavy quarks on the lattice 

Lattice QCD allows for the computation of 
physical quantities from first principle, i.e. di- 
rectly from the Lagrangian of QCD. But working 
in a finite and discrete space-time implies system- 
atic errors, which have to be under control. The 
main problem with heavy quarks is that the dis- 
cretization errors (when a traditional Wilson-like 
fermionic action is used) are proportional to pow- 
ers of the bare quark mass. Thus one has to re- 
quire am quar k < 1. To simulate a b-quark at its 
physical mass (« 5GeV) on a space-time of vo- 
lume (aiV) 4 = (2 fm) 4 , one needs N 50 points 
for each space-time dimension, which is impossi- 
ble with present computers 0. 

One can simulate various quark masses in the 
regime where the discretization errors are under 
control (say around the charm quark), and ex- 
trapolate the results to the b-quark. The prob- 
lem is that this extrapolation is done over a large 
range (because mb ~ 4m c ), and the associated 
systematic error is then difficult to control El 
Another possibility is to use effective theories, like 
heavy quark effective theory (HQET) or non rel- 
ativistic QCD (NRQCD). Since the heavy quark 
mass is much larger than the other scales (like its 
3-momentum or Aqcd), the idea is to expand the 
QCD Lagrangian in inverse powers of the heavy 
quark mass and keep only the leading terms (I 



1 Direct simulations of the c-quark are doable, but some 
care has to be taken to control the discretization errors [6] . 

2 This extrapolation can be replaced by a -more safe- in- 
terpolation, by the use of an effective theory. See [7] for a 
quenched computation of fg a done in this way. 
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give more details below). The procedure results 
an effective theory, which has to be matched with 
QCD. This matching is a source of uncertainty 
when it is done perturbatively. Recently, some 
efforts have been made to remove the largest dis- 
cretization errors by adding appropriate terms to 
the Lagrangian. One obtains a Lagrangian which 
describes both light and heavy quarks. One of the 
main problem with these methods is to compute 
the coefficients which come in front of the terms 
that one adds to the Lagrangian. 

In the next section, I present shortly these ef- 
fective theories, and their lattice discretization 
(see [8] for a pedagogical review). 

3. Effective theories 

A number of simplifications can be made when 
one deals with one or several heavy quarks. To 
make this explicit, one can derive an effective La- 
grangian, which is much simpler than the origi- 
nal. The starting point is the observation that 
the momentum of a heavy quark inside a hadron 
can be written as p — ttlqv + k. In that de- 
composition, v is the velocity, and k the resid- 
ual momentum, which is zero if the heavy quark 
is on-shell. For a heavy quark interacting with 
light degrees of freedom, the components of the 
residual momentum k are of order Aqcd, and 
therefore much smaller than m,Q. One can sepa- 
rate the higher and lower components (this refers 
to the case where the 7 matrices are written in 
the Dirac basis) of the heavy quark field Q{x) in 
h±(x) = e m Q v - x P±Q(x), where P± = ±fl Then 
one finds the effective tree-level Lagrangian[f| (see 
e.g. [9] for a simple derivation) which reads 
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where the ellipse represent higher order terms. 

The values of the coefficients given in eq. (P) 
only hold at tree level 0. They have to be renor- 



3 In principle one can also add a mass term, but it can be 
absorbed in the definition of the quark mass. 
4 In regularization scheme which preserves reparametriza- 
tion invariance, it can be shown that w kin = f for a certain 
choice of the definition of the mass 



malized to include loops effects. This is usually 
done though a perturbative matching with QCD. 
One computes a physical process at a certain or- 
der of the effective theory, and imposes the re- 
sult to be equal to its QCD value. However, on 
the lattice, a perturbative computation of these 
coefficients will lead to a results which is power 
divergent in the continuum limit. Thus a non- 
perturbative matching with QCD is needed, this 
has been done in the case of HQET [TU](I will 
give an explicit example in Sect. |4~2|) , but not for 
NRQCD. 

3.1. Heavy quark effective theory 

HQET describes a situation where only one 
heavy quark is present, like heavy- light mesons. 
It is then more natural to write everything in the 
rest frame of this heavy quark. The previous 
effective Lagrangian is seen as an expansion in 
Aqcd/tiq- With Aq CD ~ 500 MeV and toq ~ 5 
GeV the theory is expected to have roughly a 
10% precision at the leading order, and 1% at 
the next to leading order. At the leading or- 
der, when rriQ — > 00, the Lagrangian is then 
Co = h + {x) [iDo] h + (x) . It represents a heavy 
quark acting only as a static color source. One 
can see that, at this order, the light quarks are 
independent of the flavor and of the spin of the 
heavy quark. The second and the third terms ap- 
pear at the next to leading order, and represent 
respectively the interaction due to the motion and 
to the spin of the heavy quark : 
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A lattice formulation of HQET at the leading or- 
der is the so-called Eichten-Hill action [TIJ. It 
is important to note that the higher order terms 
appear only as insertion in the static green func- 
tions. Explicity, at the 1/toq order, one writes 

exp (-5iight - Shqet) = ex P ("flight - Sstat) 

X f 1 + ^ X! I^kinOkin + ^spinCspin]^ • 
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For a green function of an operator O, this means: 

(O) = (C) s tat +^ kin a 4 ^(OO ki „(^))stat 

X 

stat ? 

x 

where (0) s tat is the expectation value of O given 
by the Lagrangian Cn g ht + £stat ■ I n that way, the 
continuum limit is well defined, and the theory 
is (power-counting) renormalizable. In the pre- 
vious expressions Okin and O sp i n are the lattice 
version of the operators corresponding to © and 
([3]) respectively. The coefficients an d w sp i n 
are fixed by the matching with QCD. 

3.2. Non relativistic QCD 

In a heavy-heavy hadron, the dynamics is 
rather different than with heavy- light, and HQET 
is not an appropriate theory. Starting from 
the same effective Lagrangian ([I}, it is pos- 
sible to derive another effective theory, called 
nrqce| |12j . In a quarkonium, the typical mo- 
mentum of a heavy quark should be such that its 

-2 

kinetic and potential energy are balanced (^} ~ 
{— ). Since the distance r between the quark 
and the antiquark is of order 1/ \p\ , one find that 
(\p\) ~ a s rriQ. Then a natural separation of the 
scales is given by \v\ ~ a s : 

J) 

niQ > \p\ ~ m Q \v\ > ~ ttiqV > ... 

m Q 

Therefore in NRQCD, the power counting is dif- 
ferent from that of HQET |12ll3j . the two first 
terms of the effective Lagrangian are of order 
rriQV 2 and the third one (the 'spin' term) is of or- 
der rriQV 4 . The physical picture implies that the 
kinetic term cannot be treated as a small \jmq 
correction (as in HQET), but has to be included 
already at the leading order. On the lattice, this 
implies that the continuum limit cannot be taken. 

The authors of [12] gave an improvement pro- 
gram to reduce the errors below 10%. The prin- 
ciple is to add next to leading order corrections 
to this Lagrangian. First there are some rela- 
tivistic corrections (v 2 ~ 0.1 for the b-quark, and 

5 NRQCD can also be used for heavy light mesons 



0.3 for the c-quark). Then the leading discretiza- 
tion errors are of order 0(a 2 p 2 ) and 0(a^~^), 
which are the same order of magnitude as the 
relativistic corrections. These corrections to the 
Lagrangian introduce new coefficients, which con- 
tain power law divergences (in g 2 /amq), because 
of the non-renormalizability of the theory. They 
are computed in perturbation theory, implying 
that the lattice spacing should not be too small 
(typically a ~ 1/mg). To reduce the lattice spac- 
ing (and increase the precision), one should add 
more and more terms to the Lagrangian. Naively, 
these radiative corrections can contribute for 20 
or 40% errors, so these perturbative calculations 
with care. A lot of work has been done is the 
last years to obtain lattice results for the hadron 
spectrum with a very high accuracy, and the in- 
clusion of dynamical fermion plays a crucial role, 
as discussed in Sect. 14.11 

3.3. Relativistic heavy quarks 

The lattice formulations of both HQET and 
NRQCD are obtained from their continuum ver- 
sion. One first writes down an effective La- 
grangian in an euclidian space-time at a certain 
order, then this Lagrangian is discretized. Of 
course the resulting theory is only valid for heavy 
quarks (and low velocity for NRQCD). 

This is rather different in the Fermilab ap- 
proach (see e.g. |I4ll5lf6j ). The aim is to find a 
lattice Lagrangian able to describe both small and 
large masses. The idea is to use on-shell Symanzik 
improvement, treating both a and 1/rag as short 
distances, in such a way that the theory still 
makes sense when arriQ > 1. Starting from the 
standard Wilson-Dirac operator, one can write 

AT"" = £ W ilso„ + « dim °" 4 C00lat • (4) 

o 

The (dimension 5) kinetic and spin operators ap- 
pear in Oiat , in the the previous expression. Their 
coefficients cq are functions of arriQ which are 
tuned to reproduce their continuum values (taken 
e.g. from NRQCD). Unfortunately this match- 
ing is done done (at least partially) in perturba- 
tion theory. The claim is that when a — > 0, the 
action reduces to the (improved) Wilson action. 
Then contrary to NRQCD, the continuum limit 
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of ([I]) exists. Once all the coefficients have been 
computed at a sufficient order, the Lagrangian 
(gD, has the properties of NRQCD (or HQET) 
for aniQ > 1, and reduces to the light action 
for arriQ < 1. In the same spirit, the authors of 
[17] reduce the discretizations errors up to order 
oAqcd , by using 4 dimension 5 operators (with 4 
different coefficients). In a very recent work [18] . 
it is shown that the coefficients of these opera- 
tors are not independent, and that, indeed, all 
errors of order 0(a|p|) and 0((amq) n ) can be re- 
moved (for arbitrary n) by a proper choice of 3 
coefficients (including the bare quark mass mo). 
Moreover, contrary to NRQCD or to the Fermilab 
method, these coefficients can be computed non- 
perturbatively [19] . and then would represent a 
great improvement. 

4. Physical results 

4.1. Unquenched results 

A few years ago, in a common paper [20] . the 
MILC, HPQCD and Fermilab collaborations pre- 
sented a comparison of lattice computation of 
hadronic observables with the experimental re- 
sults. In the quenched approximation, the results 
agree within 10%, and when the effects of quark 
vacuum polarization are included, the agreement 
increases to 3%, as one can see in Fig [I] Last 
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Figure 1. Comparison of quenched and un- 
quenched results with experiments. 

year, the same collaborations achieved the com- 



putation of various quantities which were experi- 
mentally unknown or poorly known at that time, 
like the decay constant of the D s and the D + 
mesons, or the mass of the mesons. After 
precise experimental measurement of these quan- 
tities, the main conclusion is again that the inclu- 
sion of dynamical fermions is very important for 
the agreement between lattice and experimental 
results [2"T] (see e.g. Fig[2]). The b-quark was sim- 
ulated with a 0(a 2 , w 4 )-improved NRQCD action, 
and the c-quark with the Fermilab action. For the 
light dynamical fermions, the staggered action 
was used. The problem with these fermions, is 
the use of fourth root trick, to eliminate the non- 
physical 'tastes'. Since this can introduce non- 
locality, this can be potentially dangerous (for a 
recent review about this, see [22] )• The advantage 
of using staggered fermion resides in the fact that 
they are numerically very cheap. Simulating light 
quarks is then easier than than with other actions 
(in this work the sea quark masses go down to 
~ m s /6). 
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Figure 2. Comparison of quenched (nf = 0) and 
unquenched (nf = 2,3) results of fo s with its 
experimental value. 

4.2. Quenched results in HQET 

I report here some calculations in the b-sector, 
which present some theoretical advantages, de- 
spite the use of the quenched approximation. 

A computation of the bag parameter of B s — B s 
mixing in the static limit has been presented at 
this conference [53]. The main advantage of this 



5 



computation is that the light quark is simulated 
with the overlap action [24j . which presents a 
chiral symmetry at finite lattice spacing. This 
simplifies the renormalization of the 4-quark op- 
erators, and reduces the systematic errors. The 
result is fiff = 0.92(30. 

The two last exemples concern the mass of the 
b-quark and the decay constant fg a . One com- 
mon point is the use of a small volume (here of 
space extent L\ ~ 0.4 fm) to simulate a rela- 
tivistic b-quark with discretization errors under 
control. The authors of have used HQET 
together with the step scaling method, to com- 
pute these quantities [26|27j . One starts by the 
computation of a finite- volume observable Q(L\). 
The evolution to a larger volume Lj+i = kL\ is 
given by a step scaling function (ssf), defined by 
C(^2) = &(L2)/Q(Li). If the volume of space ex- 
tent Lqo = Lfinai is large enough to be considered 
as an infinite volume, one has 

^(Lo.) = <T(L &nal )...a(L 3 )a(L 2 )$(L 1 ) (5) 

In the r.h.s. of |5]), ^(Li) is computed in QCD, 
and the ssf are interpolated to the b-quark mass 
by using QCD with lower masses together with a 
static calculation. The l.h.s can be an HQET pre- 
diction or an experimental input. One choice of 
observable is the pseudo-scalar mesons mass. In 
the infinite volume, it is fixed to the experimental 
value of the B s mass. The b-quark mass is then 
extracted from the r.h.s of (5j). Of course, in the 
case of the decay constant, the l.h.s. of © is a 
prediction. The 'static' results are m^ s (mb) = 
4.421(67) GeV and F Bs = 191(6) MeV. 

A non perturbative calculation of the b-quark 
mass including the 1/niQ terms has been done 
in [55] . This is a direct application of the method 
introduced in [TU] , where a calculation of the b- 
quark mass in the static approximation is also 
performed. Here I just give a quick overview of 
the strategy. For simplicity, I start with the static 
case. At this order the meson mass is given in 



6 Here and in the following, the errors do not take into 
account the use of the quenched approximation. 



terms of the bare b-quark mass by 

m B = m ba xe + ^ (6) 

We can define an observable $(L) both in QCD 
and in HQET, such that ^(Loo) = tob- In a finite 
volume Li, the equivalent of (|6|) is 

$(L X , M b ) QCD = m bare + r stat (i!) (7) 

where r stat (ico) = E stat . Since the bare pa- 
rameters are independent of the volume, one can 
use (0 to substitute mbare m ©• This gives 
an expression where [E stat — r stat (Li)] appears. 
The terms E stat and T stat contain divergences li- 
near in the inverse lattice spacing, but they can- 
cel in the difference. However, since is much 
larger than L\, it is in practice very difficult to 
find a lattice spacing common to these volumes. 
Instead, one can introduce an intermediate vo- 
lume L2 — 2Li, and the step scaling function 

o* tat = £ 2 [r ,rt " t (L2)-r" tat (ii)] . (8) 

Then one can write: 

^.stat 

m B = M b ) QCD + [£ stat - r stat (£ 2 )l + 

Now the cancellations of the divergences occur 
in E stat - r stat (i 2 ) and in o- stat . Of course the 
procedure can be repeated, but it appears suffi- 
cient in practice to do the computation with 3 
different volumes ~ 1.4 fm > L 2 ~ 0.7 fm > 
L\ ~ 0.35 fm. Finally, the mass of the B me- 
son is fixed to its experimental value and one ex- 
tracts the RGI b-quark mass by an interpolation 
($(Li,Mb) has to be computed for a few quark 
masses around the b-quark mass). The whole 
procedure can be generalized at the next order 
of HQET. The startingpoint is to rewrite ^ in- 
cuding the 1/toq term|j 

m B — m b are + & + ^kin-^ 

Then one has to find an observable to eliminate 
cjkin, and applies the same techniques than in the 
static case (e.g. two other ssf are needed). The 

7 In principle there should be a term a; sp i n _E spln , but it 
can be eliminated by considering a spin-average B meson 
mj = ^(rriB + 3itib* )■ Then the spin-splitting becomes 
a separate issue. 
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result is m^ s (m b ) = 4.437(48) GeV. 
This non-perturbative method presents several 
good features. It is based on general hypotheses, 
and so can be applied to various quantities, like 
heavy-light decay constants. The cost is rather 
moderate, so it should be possible to generalize 
it to full QCD. Since the problem of power-law 
divergences is cured, the theory is well-defined in 
the continuum limit. Because calculations can be 
done at the NLO of HQET (the NNLO correc- 
tions are expected to be very small) , a very good 
precision can be reached. 

5. Outlook 

Heavy flavor physics on the lattice is a very ac- 
tive field, and lot of progress has been achieved 
recently. Some of them are technical (all-to-all 
propagators [29], noise reductions [30], new algo- 
rithms), and other are conceptual (new strategies, 
new effectives actions, non-perturbative renor- 
malization). Results with simulations done with 
dynamical quarks are already available. Although 
a lot of them are done with the controversial stag- 
gered quarks, simulation of dynamical fermions 
with other actions are on their way. They have 
a major role to play in constraining the standard 
model and in the search for new physics. 
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